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Abstract 

We propose TBA integral equations for multiparticle soliton (fermion) 
states in the Sine-Gordon (massive Thirring) model. This is based on T- 
system and Y-system equations, which follow from the Bethe Ansatz solution 
in the light-cone lattice formulation of the model. Even and odd charge sec- 
tors are treated on an equal footing, corresponding to periodic and twisted 
boundary conditions, respectively. The analytic properties of the Y-system 
functions are conjectured on the basis of the large volume solution of the sys- 
tem, which we find explicitly. A simple relation between the TBA Y-functions 
and the counting function variable of the alternative non-linear integral equa- 
tion (Destri-deVega equation) description of the model is given. 



1 Introduction 



A better theoretical understanding of finite size (FS) effects is one of the most 
important problems in Quantum Field Theory (QFT). The study of FS effects is a 
useful method of analysing the structure of QFT models and it is an indispensable 
tool in the numerical simulation of lattice field theories. 

Liischer derived a general formula for the FS corrections to particle masses 
in the large volume limit. This formula, which is generally applicable for any QFT 
model in any dimension, expresses the FS mass corrections in terms of an integral 
containing the forward scattering amplitude analytically continued to unphysical 
(complex) energy. It is most useful in 1 + 1 dimensional integrable models 0, where 
the scattering data are available explicitly. 

The usefulness of the study of the mass gap in finite volume is demonstrated [3j 
by the introduction of the Liischer- Weisz- Wolff running coupling that enables the 
interpolation between the large volume (non-perturbative) and the small volume 
(perturbative) regions in both two-dimensional sigma models and QCD. 

An important tool in the study of two-dimensional integrable field theories is the 
Thermodynamic Bethe Ansatz (TBA). This thermodynamical method was initiated 
by Yang and Yang ^ and allows the calculation of the free energy of the particle 
system. The calculation was applied to the XXZ model by Takahashi and Suzuki [Sj 
who derived the TBA integral equations for the free energy starting from the Bethe 
Ansatz solution of the system and using the "string hypothesis" describing the dis- 
tribution of Bethe roots. The form of the resulting TBA equations strongly depends 
on whether the anisotropy parameter po (see Section 2) is an integer, rational or 
irrational number. 

The TBA equations also determine FS effects in relativistic (Euclidean) invariant 
two-dimensional field theory models where the free energy is related to the ground 
state energy in finite volume by a modular transformation interchanging spatial 
extension and (inverse) temperature. Zamolodchikov ^ initiated the study of TBA 
equations for two-dimensional integrable models by pointing out that TBA equations 
can also be derived starting from the (dressed) Bethe Ansatz equations formulated 
directly in terms of the (infinite volume) scattering phase shifts of the particles. In 
this approach the FS dependence of the ground state energy has been studied [7] 
in many integrable models, mainly those formulated as perturbations of minimal 
conformal models. 

TBA methods have been developed also in lattice statistical physics [8^. Actually 
the basic functional relations, Y-system equations, TBA integral equations and also 
other techniques playing important role in the TBA analysis of continuum models 
have been originally introduced here. TBA equations describing the FS energy of 
excited states were proposed for some models. Also non-linear integral equations, 
similar to the Destri-deVega equations, appeared here first. An important model, the 



1 



tricritical Ising model perturbed by its 0i^3 operator, has been studied in detail p. 
The TBA equations describing all excited states are worked out in this example, 
both for the massive and the massless perturbed models. 

The TBA description of excited states is less systematic in continuum mod- 
els. The excited state TBA systems first studied [Ull Ej are not describing particle 
states, they correspond to ground states in charged sectors of the model. An interest- 
ing suggestion is to obtain excited state TBA systems by analytically continuing [TT] 
those corresponding to the ground state energy. TBA equations for scattering states 
were suggested for some perturbed field theory models by the analytic continuation 
method 12j. Excited state TBA equations were also suggested for scattering multi- 
particle states for the Sine-Gordon model at its = 2 supersymmetric point [13j. 

The BLZ program ^1] has been initiated to derive functional relations, Y- 
systems, and TBA integral equations, both for the ground state and for excited 
states, directly in the continuum. The construction is complete for the conformal 
field theory limit of the models only although some examples in massive perturbed 
models are also worked out. The original construction is based on the i?- matrix 
corresponding to the quantum deformed loop algebra Uq{A^^^) and can be applied 

(2) 

to minimal models perturbed by their 01^3 field but it works also in the Uq{A2 ) 
case pn], which is relevant to the 0i^2 and 02.1 perturbations. In this family the 
TBA equations for the first two excited states in the non- unitary perturbed A^s^s 
model are worked out explicitly. 

An alternative to the TBA equations is provided by the Destri-deVega (DdV) 
non-linear integral equations. They were originally suggested fHl for the ground 
state of the Sine-Gordon (SG) or the closely related massive Thirring (MT) model, 
but can be systematically generalized to describe excited states as well; the (com- 
mon) even charge sector of the SG/MT model [T7j, the odd charge sector in both SG 
and MT and even the states of perturbed minimal models, which can be repre- 
sented as restrictions of the SG model ^Hj- The advantage of this non-linear integral 
equation approach is that when available (SG/MT model and restrictions), it gives 
a systematic description of all excited states. The method has been generalized to 
higher rank imaginary coupling afiine Toda models |20j. 

Meanwhile an alternative derivation of the TBA equations, which is independent 
of the vahdity of the string hypothesis has been found in the lattice approach for the 
generalized Hubbard model j25 and also for the XXZ model |221- Both models are 
solvable on the lattice by the Bethe Ansatz method and the new derivation of the 
TBA equations proceeds via the introduction of T-system and Y-system functions 
satisfying functional relations which can be transformed into TBA-type integral 
equations. Instead of the string hypothesis used in the original derivation here 
the basic assumption is about the analytic structure of the Y-system functions (dis- 
tribution of their zeroes) and the validity of these assumptions can be convincingly 
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demonstrated numerically, at least for small lattices, by starting from the numer- 
ical solution of the Bethe Ansatz equations. The Takahashi- Suzuki results for the 
ground state energy are reproduced and new TBA equations for the excited states 
are found in this way. Actually for the excited state problem the TBA equations are 
supplemented by quantization conditions restricting the allowed values of particle 
momenta, which is natural for particles confined in a box. 

In this paper we borrow the ideas of the lattice approach and apply them to the 
continuum QFT case. Following the steps of the lattice construction, we systemati- 
cally build the T-system and Y-system functions for all multi-soliton (multi-fermion) 
states of the SG (MT) model. We find the appropriate continuum TBA equations 
and quantization conditions and a link between the TBA and the DdV equations. 

Our starting point is the light-cone lattice regularization of the SG/MT 
model and the Bethe Ansatz solution. We use twisted boundary conditions on 
the lattice and show that by changing the value of the parameter characterizing the 
boundary conditions we can describe both the even and the odd sectors of the model 
on an equal footing. Previously only the even sector was treated in the light-cone 
lattice formulation and the properties of the states belonging to the odd sector were 
conjectured by postulating the corresponding DdV equation in the continuum limit 
|18j . Our result for the odd sector shows that the conjectured DdV equations are 
indeed correct in this sector. 

Next we introduce the T-system and Y-system functions and write down the 
Y-system functional relations. In the even sector we recover the usual Y-system 
equations but we find that in the odd sector the Y-system equations are somewhat 
modified. Although the small modifications (with respect to the even sector) only 
occur at the end of the corresponding diagrams, this result shows that a naive 
analytic continuation of the even sector TBA equations would lead to incorrect 
results for the odd sector. 

There is a link between the TBA Y-system functions and the "counting function" 
of the DdV approach. This link is explicitly given by the simple formula (|T8|l and is 
used to write down explicitly the solution of the TBA Y-system for asymptotically 
large volume. 

To transform the Y-system equations into TBA integral equations (and quantiza- 
tion conditions) we need to know the analytic properties of the Y-system functions, 
in particular the distribution of their zeroes. From the explicit solution we know 
this distribution in the infinite volume limit. Our main assumption in this paper is 
that the qualitative properties of this distribution remain the same for finite volume. 
In the TBA language the states are characterized by specifying the source terms in 
the integral equations and since the functional form of the source terms is volume 
independent, this can be read off from the known solution in the large volume limit. 

Using the above method we write down the complete set of equations determining 
multi-particle momenta and energies in the TBA approach. We have verified the 
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validity of our main assumption by numerically comparing the results of the TBA 
approach with those obtained using the DdV equations. In a particular case (po = 4), 
where we found exact results for one-particle momenta and energies on the TBA 
side, we were able to verify the agreement by numerically solving the corresponding 
DdV equations to 20 digits. 

For simplicity, in this paper we restrict attention to the repulsive case po > 2 
of the SG coupling and consider integer pq only. Moreover, we consider multi- 
soliton states only and no states containing both solitons and anti-solitons at the 
same time. We believe that appropriate TBA systems can be found also for more 
general couplings and more general states but the special cases we are considering 
in this paper are sufficient to present the main ideas and assumptions. Finally 
we note that although in the SG/MT case the excited state TBA description is 
"superfluous" since we already have the DdV equations to study FS physics, one 
can hope that the simple pattern of the excited state TBA systems we find here 
means that similar systems can be found also for other models, where no DdV type 
alternative is available. 

The paper is organized as follows. The light-cone lattice regularization of the 
SG/MT model and the Bethe Ansatz solution is briefiy recalled in Section 2. In 
Section 3 we introduce the T-system and Y-system functions and establish the link 
between the TBA Y-system and the "counting function" of the DdV approach. 
In Section 4 simple properties of the counting function are recalled. In Section 5 
we write down the light-cone lattice TBA equations together with the quantization 
conditions. In section 6 we take the (finite volume) continuum limit of the equations. 
Multi-particle energy and momentum expressions are given in Section 7. In Section 
8 we discuss a special case (po = 4) where we can find exact expressions for the one- 
particle energies and momenta. In section 9 we discuss the DdV equation and its 
analytic continuation to the entire complex rapidity plane. In Section 10, using the 
link between the DdV equation and the TBA Y-system, we find the explicit solution 
of the infinite volume limit of the TBA problem. The ground state (no particles), 
one-particle, and two-particle problems are discussed in detail in Sections 11, 12, 
and 13, respectively. The numerical comparison of the results of the TBA approach 
with those obtained using the DdV equations is briefiy described in Section 14 
and finally our conclusions are summarized in Section 15. The technical details of 
the transformation of Y-system type functional relations into TBA type integral 
equations is discussed in the Appendix. 

2 Light-cone approach to twisted SG model 

Our starting point is the Bethe Ansatz solution of the integrable lattice regulariza- 
tion of the Sine-Gordon field theory |2S1- Here we briefiy summarize the results of 
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this approach ^1 QH| to the Sine-Gordon (and massive Thirring) model. The 
fields of the regularized theory are defined at sites ( "events" ) of a light-cone lattice 
and the dynamics of the system is defined by translations in the left and right light- 
cone directions. These are given by transfer matrices of the six-vertex model with 
anisotropy 7 and alternating inhomogeneities. This approach is particularly useful 
for calculating the finite size dependence of physical quantities. We take points 
{N even) in the spatial direction and use twisted periodic boundary conditions. The 
lattice spacing is related to L, the (dimensionful) size of the system: 



The physical states of the system are characterized by the set of Bethe roots 
{wj, j = 1, . . . ,m < which satisfy the Bethe Ansatz equations (BAE) 



Q{wj + 2i) _ To{wj + 



Q{wj - 2i) To{wj 



j = l,...,m. 



Here 



N/2 



where G is the inhomogeneity parameter, 

m = 

and 



sinh 1^ 
sin 7 



(2) 



(3) 



(4) 



(5) 



where the parameter u is characterizing the twist of the boundary condition. We 
shall also use the parametrization 
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TT 



Po 



IT 



P+l 



(6) 



Finally the energy (E) and momentum (P) of the physical state can be obtained 
from the eigenvalues of the lightcone transfer matrices: 



Ja{E±P) 



+2iuj 



(7) 



Besides the usual procedure, taking the thermodynamic limit (A^ —>■ 00) first, 
followed by the continuum limit (a — > 0) one can also study continuum limit in finite 
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volume by taking N ^ oo and tuning the inhomogeneity parameter B simultane- 
ously as 

2 2N 

e = ln-f = ln^, (8) 
Ma I 

where the mass parameter M is the infinite volume physical mass of the Sine-Gordon 
solitons and we have introduced / = ML, the dimensionless size of the system. 



3 T-system and Y-system 

Besides Q{^) and To{^) it is useful to introduce the (half-)infinite sequence of 
functions Ts(^) for s = — 1, 0, 1, . . . In the six-vertex model these are eigenvalues of 
the transfer matrices corresponding to higher spin representations (in the auxiliary 
space) and are obtained by the fusion procedure. Equivalently we can simply define 
them as T_i(^) = and 

s 

TsiO = Q{^ + ^ + ^s)Q{^ -t-ts)Y,^[^ + ^(2J - ^)] (9) 

j=0 

for s = 0, 1, . . . , where 

It is obvious from the definition that the Tg are periodic functions, 

Ts{^ + 2tpo)=Ts{0. (11) 

It can also be verified that they satisfy the T-system equations [22] 

T,{^ + t)Ts{^-t)=Ts+i{OTs-i{0+To{^ + ^ + ^s)To{^-t-ts), s = 0, 1, . . . (12) 

which can then be used to recursively calculate them starting from Ti(^). 
The next step is to introduce the functions 

v(n- ^s+i(Or,_i(0 

or equivalently 

As consequence of the T-system equations (fT^ they satisfy the Y-system equations 
+ -^) = [l + >^.+i(0][l + ^^.-1(0], s = l,2,... (15) 
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Finally we introduce the exponential counting function 



Tie; + 2-0 ro(? - -0 ,,„, 
^«' = QK-2.)r„K + .)- <^«> 

which can be used to rewrite the Bethe Ansatz equations as 

^K) = -l j = l,...,m. (17) 

Combining the above definitions it is easy to see that the Y-functions can simply be 
expressed in terms of the exponential counting function starting with 



i + ri(O = [i + -^(e + 0] 



1 

1 + 



^(e - ^) 



and then using the Y-system equations (fTH|l recursively. 

The above algebraic considerations are supplemented by the observation that 
although the T!,(0 as defined by Q and (fTUj) appear to be trigonometric rational 
functions but in fact they are trigonometric polynomials. This is a consequence of 
the Bethe Ansatz equations (j2|). It is also easy to see that 

TsiO ~ ei^l«l (19) 

asymptotically for C, —>■ ±oo. The Y-functions are trigonometric rational functions, 
bounded for large |,^|. 

In the rest of this paper we restrict our attention to the following special case of 
the problem. 

• Po = P + 1 integer (p > 2) 

• Wj all real and different 

• UJ = OT UJ = ^ 

We have chosen these restrictions because while they make the discussion technically 
much easier than for the general case they still contain many physically interesting 
cases. We will discuss in detail 

• case A: a; = 0, 

• case Si: = |, po odd, 

• case B2: uj = ^, pq even. 
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For the special case po integer the T-functions satisfy 

T,M-Tpo-2{0 = 2{-ircos{u;po)Toi^ + tpo) (20) 
and for cases A and B2 this allows us to define 

A-K) = (-ircos(.p„)5|^. (21) 

which, together with {Fs(^), s = 1,2, . . . ,po — 2} form a finite Y-system of type 
Dpg, because 

Ki^ + t)Ki^-t) = 1 + 1^0-2(0, (22) 
[l + K{Of = l + i;o-i(0- (23) 

For case Bi implies the reflection symmetry 

Tpo-i-kiO = rpo-m(0 k = 0,l,...,po (24) 
Yp,-i-k{0 = Ypo-i+kiO k = 0,l,...,po-l (25) 

and this symmetry allows us to truncate the Y-system after i^o_i. We call the 
corresponding finite Y-system of type A2pg_^. 



1 2 p-2 p-l P 



Figure 1: Dynkin-diagrams associated with Dp+i and ^2p-i type Y-systems 



To analyse the Y-system equations further we need some definitions. Let us 
write 

N = 2{m + d), (26) 

where < d < y- A special case of A, for which also d < po will play an important 
role in the following. We call this sub-case A. Similarly we denote by B2 the sub-case 
of B2 with d<f.We next define C(0 by 

po-l 

C(0 = E^[^ + ^(2j + l-po)]. (27) 
3=0 



8 



It is easy to see that for cases A and B2 

C(e + 20 = C(0 (28) 

and 

Tpo-iiO = C(OQ(e + iPo)Q{^ - iPo)- (29) 

Further 

TpoiO = i-ly''KToi^ + ^Po)+a^~^)Qi^ + ^ + ^Po)Q{^-^-^Po), (30) 
rp,_2(0 = {-ir+'nToi^ + ipo) + C(e - ^)Q{^ + I + ^Po)Q{^ - I - ^Po), (31) 

where k = cos{ujpo). 

• — o — o o — o 

1 2 3 p-3 p-2 

Figure 2: Dynkin-diagram associated with Ap_2 type Y-systems. 

Because C(0 periodic with period 2i and is free of singularities in the strip 
|Im ^1 < 1 as can be seen from (j^^ . it is actually an entire function and thus a 
Laurent polynomial in the variable e'^^. U d < pq, as is the case for A, then the only 
such polynomial consistent with the asymptotic equation (jTHl) is C(0 = co = const. 
Similarly for B2 ()19j) forces ({C,) to vanish. In this latter case the Dp^^ type Y-system 
is further reduced since from we have in this case 

T,„_i(O = and Yp„„2(O=0. (32) 

The remaining Y-system for {Fs(^), s = 1, 2, . . . ,po — 3} is of type Ap^^^. 
For later use we note that in case A 

Tpo_i(0 = coQ(e + tPo)Qi^ - tpo) (33) 

and this implies that Tp(,„i(^) has no zeroes in the strip |Im ^| < 1. Similarly in 
case B2 

7^0-2(0 = (-l)™+^(-l)^To(e + tpo) (34) 

implies that Tp^^_2{^) has no zeroes in the strip |Im ^| < 1. 

To summarize, for integer po = p + 1 we found three different types of finite Y- 
systems. The Y-functions can be associated with the nodes of the Dynkin-diagrams 
shown in Figures 1 and 2. Cases A and B2 correspond to the Dp+i-type diagram of 
Figure la and case Bi corresponds to the A|p_]^-type diagram of Figure 16. The Ap_2- 
type diagram of Figure 2 corresponds to the sub-case B2. The Y-system equations 
are of the form 

Wai^ + t)Wai^ - = n [1 + ^biOY"' ■ (35) 

6 
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Here Wa{C,) = YaiO almost all cases with the exception of the last two nodes of 
the -Dp+i diagram (Figure la). In this case 

W,+,iO = W,{0 = KiO- (36) 

The matrix elements lab are defined as follows. lab is zero if nodes a and b are not 
connected by links and it is unity if the nodes are connected by a single line. Finally 
the oriented double line at the end of the A|p_i type diagram means 

Ip-i p = 1 p-i = 2. (37) 

Our -Dp+i and y4p_2-type Y-system equations coincide with the usual Y-system 
equations [21] associated with the Dynkin-diagrams of these simply laced Lie al- 
gebras. Our equations corresponding to the diagram of Figure 16 correspond to 
those solutions of the standard Y-system equations for the Lie algebra A2p-i that 
are symmetric under reflection of the Dynkin diagram. 



4 The counting function 

The next step in the Bethe Ansatz solution of the model is the definition of the 
"counting" function Zjsf [ini HZl EE] , which will play an important role in our consid- 
erations. In we have already defined the exponential counting function which 
can be used to find the positions of Bethe Ansatz "holes", real solutions of (fT7|) . 
different from the Bethe roots: 



J^(^„) = -1 a = l,...,H. (38) 

To define Zjsr properly, we need some more definitions. If the function f{z) is 
analytic and non-vanishing (ANZ) in some simply connected domain U then we 
can define its "logarithm" log/(2;) as the primitive function of f'{z)/f{z), up to a 
constant, which can be fixed by giving the value of log f{zo) for some zq G U. 

Using the above definition we now introduce the functions 

sinh — {xT — ^1 

0,(0 = log —-4)-—^, 0.(0) = (39) 

smh ^ («r + 4 ) 

in the strip jlm^ < t (for < r < po). For later use we note that 

0,(±cx)) = ±7r^^. (40) 
Po 
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The counting function Z^- is now defined as^ 



TT 



m 



M^-w,). (41) 



The precise relation between Z^r and the exponential counting function is 

^(^) = (-l)^e'^^(^«), (42) 

where 

6 = m+— (mod 2). (43) 



TT 



Using pUj) we can calculate 



Zn(±oo) = ±n {m + 2d- — ] . 

PoJ 



(44) 



In the following we will assume that the Bethe Ansatz solution we are considering 
is such that the counting function is essentially monotonically increasing. We require 
that 



Z' 



N \ 77^ 



> 

> 

> 



J = i, . . . ,m, 
a = 1, . . . , H, 
for large |^| 



(45) 
(46) 
(47) 



and call such solutions non-degenerate. We will see that all solutions are non- 
degenerate for large enough I but we assume that this property of the solutions is 
sustained also for smaller values of I. 

For non-degenerate solutions (jl^ can be used to obtain a simple formula that 
gives the number of holes in terms of Bethe Ansatz data: 



CO 



UJ 



-- : 

71 

2 ■ 



H 
H 



2d -2 
2d - 1-2 



1 d 



d_ 

Po 



5 = m(mod2), (48) 

6 = m + l (mod 2). (49) 



(In the above formulae the square brackets stand for integer part.) 

(jIHj) and (Uni) show that the u = sector originally studied [121 HI] corresponds 
to states containing an even number of holes, whereas states with odd number 

"'^The ^ factor is present in the argument of Zjv because our variable ^ differs by this factor 
from the usual rapidity variable. 
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of holes are in the uj = ^ sector. The twist parameter u has been used in the 
description of restricted Sine-Gordon models [12] but here we see that it can also 
be used to study the odd particle sector in the original (unrestricted) Sine-Gordon 
and Massive Thirring models. The description of this important sector of the theory 
(containing the physically most interesting one-particle states) has been based on 
conjectures [TH] but as we see here it can be studied on an equal footing with the 
even sector. 



5 Lattice TBA equations 

We now translate the the Y-system functional relations ()35p into TBA integral 
equations As explained in the Appendix, we need to know the zeroes and 

poles of Wa{^) in the strip |Im^| < 1 together with their asymptotic behaviour. We 
will call this strip the main strip. From the definitions (|T!I|) and 1)2111 it is clear that 
all Wa{^) are bounded at infinity and that in the main strip they can have zeroes 
only, there are no poles. It is also clear from these definitions that the zeroes of 
Wa{^) are inherited from the neighbouring Ts{^) functions. 
Let the set of zeroes of Tk{^) in the main strip be 



h = {yP} " . (50) 

I. J n=l 

The zeroes are not necessarily different from each other, for example 

n..„. (51) 



2 y n Ti ^ 
The zeroes of Ti(^) are solutions of the equation 

Hy^^) = -1, (52) 

different from the Bethe roots Wj. This set includes the holes defined in (jHH|) and 
may contain C additional complex "holes", complex solutions of ()52|) in the main 
strip. Thus, Ri = H + C. 

The set of zeroes of Wa{^) will be denoted by 

ga = {e)}f:, (53) 
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and is given by 



Dp+i : Qa 

% 

Ap-2 : (ja 

Qp~2 



fa- 


-1 




-1 


Qp^ 


-1 


= fp. 


-1 


fa- 


1 




-1 


fp- 


'1 


Ufp4 


-1 = 


Ta- 


-1 




-1 


fp~ 


-3 







1 (54) 
(55) 

1 (56) 
(57) 

3 (58) 
(59) 



We note that in we have used the fact that -Rp_i = as estabhshed after 
and that in the special case A Rp = Q a.s follows from (jHSj)- 

To be able to use the results of the Appendix we make the following definitions. 

Uu) = Y,IabH^ + Wi,{u)l (60) 
h 

aM-lvjy-Jg-^y (62) 

(In the last formula V stands for principal value integration.) 

Using the main result ()17()|1 of the Appendix we now transform the Y-system 
equations into TBA integral equations: 

Qa 

WaiO = ^allr{^- 4")) exp {/5,(0} , (63) 

where da is the sign of Wa{oo) and r(^) is defined in fllTlj) . 

The integral equations are not sufficient to reconstruct the functions Wa{^) 
completely. They have to be supplemented by the "quantization conditions" 

1 + W^a(yi"^±0 =0 n = l,...,Ra a = l,...,p. (QA) 

(a = 1, ... ,p — 2 for the Ap_2 system.) 

It is easy to see that (jMj) follows from the Y-system equations but it is also true 
that (j63|l and (j64|) together are completely equivalent to the Y-system equations 
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6 Continuum limit 

We are interested in the finite volume continuum limit, which is obtained by taking 
N ^ oo and tuning the inhomogeneity parameter simultaneously according to 
In the Bethe Ansatz language in this limit we have to take m — >■ cxd while keeping 
d, H fixed. 

We assume that the Y-system functions Wa{^) have well-defined continuum lim- 
its. The corresponding limit of the set of zeroes will be denoted by r^: 



h = {yr\ , — r, = W:'\ , (65) 

k J n=l K J n=l 



where Rk Rk, since some of the zeroes can disappear in the process by going to 
infinity. 

Actually, all the zeroes of To(^) go away in the limit and the corresponding factor 



N_ 
2 



c„(0^n.(«-,r)^(-i)*{}^}^| (aa) 

becomes 

(_l)f e-'™^^5« (67) 

in the continuum limit. 

We denote by cTq the sign of Wa{oo) in the continuum limit, which can also differ 
from (Tq. 

In the continuum limit 

Qa ^ 1a={zi-^}t,^ (68) 

where 

qi = r2, qa = ra-iUra+i a = 2,...,p-l (69) 

(a = 2, . . . , p — 3 for the Ap_2 case) and 

1p = Qp+i = ^p-i (-Dp+i), 

q, = 2-rp_i iA%_,), (70) 
qp-2 = rp_s {Ap^2 P > 5). 

With these definitions the TBA integral equations in the continuum limit become 

Qa 

W^a(0 =^ae-^-'™^^5« J] r (C-zi")) exp{/3.(0}, (71) 

a=l 
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which have to be supplemented by the quantization conditions 

l + W^a(z/i"^±^) =0 n = l,...,Ra a=l,...,p. (72) 

(a = 1, ... ,p — 2 for the v4p_2 system.) The exponential factor e"'™^*^!^ is present in 
the TBA equation for a = 1 only. This is indicated in Figures 1 and 2 by colouring 
the corresponding nodes black. 

An important special case is when all zeroes are real. In this case the modulus 
of Wa (ya^'* ±i) is automatically equal to unity and (ff^ can be rewritten as 



exp -^ j^.i/sinh - a, (yj^^) + (l/i"^ - | = (73) 

for n = 1, . . . ,Ra and a = 1, . . . ,p (a = 1, . . . , p — 2 for v4p_2) . In (f73|) the notation 

7(m) = 2 arctan (r(u)) (74) 
is used. Note that |7(m)| < | for real u. 

7 Energy and momentum 

The formulae ((Zj) determinig the energy and momentum of the Bethe Ansatz state 
can be rewritten in terms of To(.^) and Ti(^) as 

T, (20 _ i) Ti (i - ^) Ti (^ - i) 

-2iaE ^ TT ''j j -2iaP ^ -Aiuj M tt '' > /yp-N 

T2(f-2.) ' Ti(z-f)- ^^'^ 

From (fT^ we see that T'i(^) satisfies 

Ti(^ + z)Ti(^ - = [1 + ^i(0] ^o(e + 20To(^ - 20, (76) 



which is of the form ()165p and can be solved by the techniques explained in the 
Appendix. For this purpose we first define 

Ti(0 = afi(0<f^/^(0, (77) 

where a is a sign, Ti(^) has the same zeroes as Ti(^), has no poles, it is positive for 
^ — > cx), and satisfies 

f,{i + i)f,{i-i) = l + Yr{i). (78) 
In (I77jl the factor $(^)^/2 solves the To-dependent part of (fTHjl : 

*({)^exp{A(0}. A(0^j/j„-J^. (79) 
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where 



B{u, \) = k{u + \ + 2i) k{u + \- 2i) k{u - \ + 2i) k{u - \ - 2i) 
Now we take the logarithm of (f73j) and get 



E 
P 



Eq + E 
Po + P 



E = S+ + S. 
P = S,-S. 



where 



2a \ 71 



S- = —ln 
2a 



[-iff^U- — 



29 



TT 



and the energy and momentum constants are 



Po 



-No 
a 



uj iN 



a 



2a 



- + ^ A — 
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«vr , . , . , /2e 

h In sm 7 — m smh 7 — 

2 y TT 



- + -m. 

a a 



^0) 



(81) 
(821 



(83) 

,(84) 
(85) 



Here the choice of the integers A^o and A^i (which are in principle arbitrary) is part 
of the regularization scheme. If we choose A^^o = A^^i = we get 



P - 

-TO — — 

a 



and 



^0 = + 7^X0 — 

a 2a \ TT 



(86) 



The precise definition of the function XoiO be given later. What is important 
here is that Eq and Pq are universal in the sense that they are the same for all 
states in a given sector characterized hj uj = or uj = ^. Omitting these unphysical 
constants leaves us with the physical energy and momentum eigenvalues E and P. 
Using the method outlined in the Appendix we can solve ()78p : 



m) = n - n <^ - ^/3) exp - / du 

a=l 0=1 ^ 



ln[l + yi(M)] 
cosh f (^ — u) 



^7) 



Here flp are complex solutions of ()52|1 . Using the definitions (jH!I|) finally we get in 
the continuum limit 



E 



M 



M 



cosh — — h cosh ■ 

a=l (3=1 
' H , C 



sinh — - + sinh 



2 



1 f°° TTU 

- du cosh — ln[l + yi(M)] 
4 J -00 ^ 



0=1 



f3=l 



TTU 



du sinh — ln[l + ^[(m)] 



^9) 
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8 The Ai case 



One of the simplest cases is p = 3, d = if = 1. In this case we can determine the 
energy and momentum eigenvalues exactly. 

In this case 5 = f and from (0 it follows that Yi(oo) > 0. Further = 

and i?2 = here and (j63|) reduces to 

Ym = W,{i)=CM{i). (90) 

It is easy to see that there cannot be any complex holes in this case. Thus Ri = 1 
and we introduce = h for the position of the single (real) hole. The solution of 
(IMjl gives 

h = h + ^^+0{N-^), (91) 

where h is the solution of 

I sinh -h = 27rMi, Mi = (mod 1) (92) 

and 

sinh (^irh) , ^ 

h2 = - ^ f . (93) 

OTT cosh (^vr/ij 

In the continuum limit we get 

Y,{0 = i-'^Y e-^""^^^ (94) 
and the energy and momentum eigenvalues are 

E = A^cosh^-^ rfMcosh^ln[l + (-l)^e-'^°^'^5-], (95) 

_ , tt/i 27rMi , , 

P = Msinh — = (96) 
Z Ij 

It is interesting to note that the momentum eigenvalues (j96|) are exactly the same 
as they would be in a free theory. The p = 3 model is not free however since the SG 
S-matrix is non-trivial at jj"^ = Qn. (The free fermion point is = 47r or p = 1.) The 
energy eigenvalues differ from those of the free theory by the second term in (jHS))- 
On the other hand, this shift is independent of the momentum quantum number 
Ml. We note that it is only the one-particle states that can be exactly solved for 
p = 3. We also note that it is not obvious in the DdV language that p = 3, H = 1 
is a special case and we were able to solve the corresponding DdV equations only 
numerically. 
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9 Destri-deVega equation 



In the Sine-Gordon (MT) model a non-linear integral equation (Destri-deVega equa- 
tion) is used to describe the Bethe Ansatz states and calculate their energy and 
momentum El CB]- This is formulated in terms of the counting function Z^- 
and its exponential, ^(0- 

Let us chose the parameter B so that there are no complex holes in the strip 
< Im^ < B. Then we can define 



/:+(0 = iog [1 + ^(0] 

We define further the function xiO by 
where 



< Im^ < B. 



X(0) = 0, 
sinh(p — l)k 



(97) 

(98) 
(99) 



2 cosh k sinh pk 

The function Xo(0 occurring in (|HH|l is defined analogously with p replaced by 

Po = P + 1- 

For non-degenerate states 



£+(0 = L+(0 = In 
and the DdV equation can be written as 



1 + 



(100) 



7r 



Zm ( tt'C ) = — ^ arctan 



arctan 



H 

sinh (l^ + e)] } + 5^x(^-/^c 



a=l 



1 r°° 

+ / du\G(^-u-ir])LJu + iri) -G(^-u + iri)L*Ju + ir])], (101) 

where < r] < B. Analogously to the quantization conditions (j6l|l supplementing 
the TBA integral equations the DdV equations are supplemented by the quantization 
conditions 



Zn { ^hr 



27rM„ 



M„ = (mod 11 



a = 1, 



H. (102) 



The counting function Z^v has continuum limit Z in terms of which the continuum 
DdV equation reads 

H 

Zge) =/sinh|e + E^(^-^") 



1 [°° 

H / (i-u [(^(^ — M — i77)Z/+('u + i?7) — (^(^ — M + 2?7)L*_(m + i?])] . 



(103) 
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where 

= In [l + . (104) 

The energy and momentum in the continuum can be expressed in terms of the 
positions of the holes ha and L+: 

E = A^^^cosh— — ^ TT^^ / du sinh —{u + irj) L^{u + irj), (105) 

P = TW^^sinh— — ^ TT^^ / du cosh —{u + ir]) L+{u + irj). (106) 



For large physical size / the function L+ is exponentially small and the integrals 
containing it in the DdV equation ()103|1 and also in the energy and momentum 
expressions can be neglected. In this approximation the quantization conditions 
()102|) are nothing but the (dressed) Bethe Ansatz equations for solitons/fermions in 
the SG/MT model. This leads to the identification 

SG : (soliton - soliton) S{e) = - exp j^x ( ^ ) [ (107) 



and 



5 = H (mod 2) 



MT : (fermion — fermion) S(6) = exp {ixi — ) \ (10^ 



TT 



5 = 0, 



where S{9) is the scattering phase for soliton/fermion scattering. 
The parameter p is related to the usual SG coupling f3 by 

p = 1 corresponds to the free fermion point and p = oo to the XY model (asymp- 
totically free point). 

For large / the first term on the right hand side of ()103|) dominates and this shows 
that Z is monotonically increasing for real ^. We already used the assumption 
that this property remains valid also for smaller /. 

The DdV equation (jlOHj) is suitable for numerical calculations. To study the 
analytic properties of the solution it is useful to rewrite it in the form 

Z = I sinh l^ + Yl -ha)+ du Gi^ - u)Q{u), (110) 
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where ^ 

Q{u) = — lim lmL^{u + irj). (HI) 

The representation ()110j) is vahd in the strip |Im^| < 2. Note that we are considering 
only the repulsive case p > 1 in this paper. In the attractive regime p < 1 some of 
our expressions have to be modified. 

We now analytically continue the counting function to the whole complex ^ plane. 
For this purpose we define 

7i(0 = -^log^i;(0 liitp) = vr l<lm^<2p-l, (112) 

72(0 = -^log^(0 72{-ip) = vr -2p + l < Im^ < -1, (113) 

where 

sinhf (^ -n 

w{0 = — r^^T TV- (114) 

"-^^ sinh|;(z + ^ ^ 

The analytical continuation of the counting function in the strip 2 < Im,^ < 2p 

is 

poo 

ZA-A=S7c + J27ii^-^-ha)+ / duK,{^-u)Q{u), (115) 

^ a=l 

where 

7^1(0 = ^(0 + ^(^-2^), (116) 

which is analytic in this strip. 

Similarly in the strip —2p < Im^ < —2 the analytical continuation is 

^ poo 

Z2(-q=Sn + J2l2{^ + i-ha)+ / duK2i^-u)Qiu), (117) 

^ a=l 

where 

^2(0 = G(0 + G(e + 22), (118) 

which is analytic in this strip. 

The exponential counting function is thus defined as 

{(_l)5gi2i(f?) 2<Im^<2p 
(_l)^e*^(50 |ImO<2 (119) 

^_lY^iZ2{^i) -2p<Im^<-2 

and it is easy to see that this defines a meromorphic function with properties 

+ 2^Po) = ^(0, = j^y (120) 

It has poles {00 many) at ^ = wj + 2i and zeroes (00 many) at = wj — 2i. 
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10 Infinite volume solution 

We have seen that for large / the integrals give exponentially small contribution to 
the counting function Z . This is also true for Z\ and Z2 although this cannot be 
seen immediately from ()115|) and ()117|) . However, it is possible to write Zi and Z^ 
in a form analogous to ()103p and show that the corresponding integrals are indeed 
exponentially small. Neglecting the integrals we have 

= [1 + C (e-^')] • <^ e*'^'°^^i«M(0 |Im^| < 2 (121) 

\ u{i)u{i + 2i) -2p < Im^ < -2 

where /C is Im ^-dependent and positive. The function m(^) is defined by 

H 

a=l 

where 

(t(^) = e*^(^) (123) 
is meromorphic with poles and zeroes on the imaginary axis and satisfies 

a{^ + tM^-i)=w{0. (124) 

Since we can express the Y-system functions in terms of the counting function 
starting with ()18|1 we are able to calculate also the Y-functions with exponential 
precision for large /. For li(^) we find 

( ri2{i-i) 3 < Im^ < 2j9 - 1 

Yi{i) =[l + {e-'^^)] ■ \ A(Oe-''=°'^^« |Im^| < 3 (125) 

\ri2{i + i) -2p + l<Im^<-3 

and using the Y-system equations (fT3jl recursively we get for = 2, . . . ,p 

n(0 =[l + (e-'^')] r7fc(0 llmei < k. (126) 

In dng) 

X[i)=u{i + i) + -^ (127) 

and 

r]2{^) = \{i + i)\{i-i)-l. (128) 
The functions rik{i) satisfy the Y-system equations 

+ i)Vk{^ - = [1 + %+i(0] [1 + %-i(0] (129) 
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for k = 2,3, . . . and this determines^ VkiO for k > 2. 

It is possible to find tlie solution of ()129p explicitly. We note first that there is 
a class of solutions depending on a parameter q and a function -8(0- Using these 
input data we first define 

fc-i 

to(0 = 0, t,(0 = Xl^'^[^ + ^(^-l-2j)] k = l,2... (130) 

j=0 

and then it is easy to show that 

Vk{^}-q + fc-1,2,... (131) 

solve ()129|) . (jl3ip is quite analogous to (fT^ and it is also true that 



The actual solution entering p25|) and ()126|) corresponds to the choice 

H 

I sinh 

2p 



q=i-lf 5(0 = n sinh :^(e-M. (133) 



This can be verified by computing ri2{^) from p28|) and from ()13ip with the above 
choice and showing that they coincide. 

Similarly to the general Y-system equations, the set of ?7-functions can also be 
truncated. In cases A and B2 we can define 

and then the relations 

<^ + -0 = 1 + vp-iiO 1 + VpiO = [1 + <0f (135) 

can be used to show that the truncated system is of type Dp^i. U H < pis also satis- 
fied {B2 case) then tp{^) = implies ?7p_i(0 = and this corresponds to the further 
truncation to Ap_2. Finally in the Bi case reflection relations analogous to (j^ 
are satisfied by tfc(0 ^^d ?7fc(0 ^^d the truncation is of type ^2^-1 accordingly. 

^ViiO = by definition. 
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Similarly to the full Y-system where the zeroes are determined by those of the 
T-system, the zeroes of tfc(^) determine the zeroes of i]k{^)- More concretely, 



A(0 : ANZ |Im^| < 3 
r]kiO ■■ ANZ |Im^| < k 
k(0 : ANZ |Im^|<p 

For later use we note that 

A(±cx)) = 2(-l)'^cos 
We now consider some simple examples. 



except for zeroes of ^2(0 5 
except for zeroes of tfc+i(0^fc-i(05 (136) 
except for zeroes of 



(137) 







[■" 





H = 



In this case we have a -Dp+i system for p >2. From ()122|) and ()127|) we get 

A(0=2(-l)^ (138) 
where, according to ()1()7|1 and p()8|l only the choice 5 = is physical. Further 

5(0 = 1, tk{0 = k, Vk{0 = k'-l, k{0=P-1. (139) 



H 



Here, according to ()107pi08|l we have two choices: 

6 = 1 for SG, 5 = for MT (140) 

and we have an A2p_i system for p > 2 even and an Ap^2 system for p > 3 odd. 



The position of the hole is denoted by hi = h. From 
case 



and p27|) we find in this 



m 



and from 



we have 



B{0 



71 



sinh — — h) 



^ ^ cosh ^ 



k odd 
k even 



:i4ii 



(142) 



(143) 



H 



Here we have a -Dp+i system for p > 2 again and only the choice 5 = is physical. 
For simplicity we restrict our attention to the symmetric case 



hi = h 



ho 



(144) 
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which corresponds to 



5^ = -cosh— - -cosh — . 145 



Now we find 

A(^) = ^Q^x{i-^i+i)+ix{£.+h+i) ^-ix(^-h-i)-ix{i+h-i)^^ (146) 

and 

1 sin (y) ^ 1^ 

tkiO = o cosh ^ - - cosh — . 147 

2 sin U j P 2 p 

11 Ground state 

The ground state of the system {H = 0) belongs to the A case and corresponds to 
a -Dp+i type Y-system for p > 2. The only physical choice is 5 = 0. From ()139p we 
see that for / — > oo there are no zeroes and all Wa{^) functions are positive. Our 
assumption is that these qualitative properties of the solution remain valid also for 
finite / values. If this is true then ()71|1 becomes 

WaiO = e-^-'^^^^^^f « exp {/3aiO} a = 1, . . . ,p + 1 (148) 

and there are no quantization conditions in this case. By taking the logarithm 
on both sides of ()148|) we reproduce the customary TBA system describing the 
ground state of the SG/MT model. Noting that both Wa{^) and /9a(0 ^-re now even 
functions we find from ()88|1 and (jHn|) 

P(°) = 0, = r du cosh ^ In [1 + Wi{u)\ . (149) 

4 J-oo ^ 



12 One-particle states 

For one-particle states both 6 = 1 (SG model) and 5 = (MT model) are possible. 
Moreover, depending on the parity of p, we have an Y-system {p > 2 even) or 

Ap-2 Y-system (p > 5 odd^). Here we only consider the SG model with /ii = in 
detail, because in this case the quantization conditions are satisfied automatically, 
which simplifies the problem enormously. 

Our main assumption will be again that the qualitative features of the / oo 
solution are also valid for finite /. From ()14H1 and p4H|l we see that 

■^We have already discussed the Ai case. 
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• W^a(0 • zeroes for a odd 

double zero at ,^ = for a even 

• all Wa{C) are negative, except Wp{^) in the v4|p_^ case 
and this leads to the TBA equations 

Wa{i) = (-1)1+^- ^-Sarlcos^^,^ [r(0]'+^"'^'' exp {[3am « = 1, . . . ,J9. (150) 

(a = 1, . . . ,]9 — 2 for Ap_2.) /?a(0 even functions and Q;a(0 ^^re odd. 

It follows that the quantization conditions (fTHj) are automatically satisfied. 
Finally from (jHHI) and (jEHl) we get 



P(^) = 0, =M-— [ du cosh — In [1 + . 

4 2 ^ 



(151) 



13 Two-particle states 

The two-particle states are the only ones among our examples where the quantization 
conditions (f7S|) give nontrivial constraints. The H = 2 states belong to the A case 
and correspond to a -Dp+i type Y-system (p > 2) and 5 = 0. 

Our conjecture is based again on the large / solution. For simplicity we consider 
only the following symmetric distribution of zeroes: 

rk = {Hk,-Hk} k = l,...,p-l, (152) 

where by convention Hk > 0, the positions of the holes are hi = Hi and /12 = —Hi 
and the signs at infinity are 

(Ti = -1, aa = l a = 2,...,p-2, crp_i = dp = o-p+i = -1. (153) 

It is important to note that ()153|1 and all subsequent formulae are valid for p > 3 
only. The important special case p = 2 will be discussed at the end of this section 
separately. 

The TBA integral equations are 

WiiO = -e-'™^^f«r(^-iJ2)r(^ + iJ2)exp{A(0}, (154) 

WaiO = r(e-/J„_i)^(e + ^a-i)r(e-i^a+i)r(e + i/a+i)exp{/3,(0},(155) 

W^p_i(0 = -r(e-ifp-2)r(e + //p-2)exp{/5p„i(0}, (156) 

WpiO = Wp+i{0 = -r{^-Hp-i)T{^ + Hp^i)exp{Pp{0}- (157) 



25 



In ()155p a = 2, ... ,p — 2 and because of the symmetric distribution of zeroes the 
functions H^a(0 PaiO even and aa{C) are odd. For the same reason it is 
sufficient to require that (|T3|l are satisfied for the positive zeroes: 

/sinh^ifi+7(/Ji-if2)+7(^i + ^2)-«i(i^i) = 27rMi, (158) 

7(ifa-i^a-l)+7(^a + i^a-l) 

+ ^(Ha-Ha+i)+iiHa + Ha+i)-aaiHa) = 27rM„ (159) 
7(iJp_i-ifp_2)+7(^P-i + ^P-2)-«p-i(i/p-i) = 27rMp_i. (160) 

In fll59p a = 2, . . . ,p — 2 and all the quantum numbers Ma are half-integers. The 
value of Ml is part of the input data, whereas all others {Ma for a = 2, . . . ,p — 1) 
are determined from consistency. Numerically we found in all cases we considered 
that Ma = 1/2 for a = 2, . . . ,p — 1. Finally the two-particle momentum and energy 
are 

p(2) _ ^(2) _ 2M cosh -Hi - — / du cosh — In [1 + Wi{u)] . (161) 

We end this section by the discussion of the special case p = 2 corresponding to 
the supersymmetric point of the SG model. Here the TEA equations are 

WiiO = e-'™^^i«exp{/5i(0}, (162) 
W^iO = Ws{0 = -r{^-HiM^ + Hi)ew{m)} (163) 

and there is just one quantization condition: 

/sinh -Hi - ai{Hi) = 2ttMi. (164) 
2 



14 Numerical iteration of the TBA equations 

We have calculated the energy and momentum eigenvalues of the ground state and 
also of one- and two-particle states from the TBA equations numerically and com- 
pared the results with those obtained using the DdV equations. We have found 
agreement to many digits in all cases we considered and concluded that our TBA 
equations are completely equivalent to the DdV equations. 

Our numerical checks covered a rather wide range of the physical size / (from 
0{1) down to / ~ 10~^) and we considered a number of different choices for the 
parameter p, including the limit p = oo. The agreement between TBA and DdV is 
especially spectacular for the Ai case where we have exact results on the TBA side. 

The infinite volume solution we found in Section 10, which we have used to help 
postulating qualitative properties of the Y-system elements, serve also as starting 
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point of the numerical iteration procedure. We found that the convergence of the 
numerical iteration is rather fast and already the starting point is a surprisingly 
good approximation not only for large volumes but also for smaller I values. In the 
H = 2 case this is true in particular for the positions of zeroes, which change very 
little during the iteration. 

15 Summary 

In this paper we proposed TBA systems describing multi-soliton states of the Sine- 
Gordon (massive Thirring) model. We derived the TBA equations starting from the 
Bethe Ansatz solution of the light-cone lattice formulation of the model, which allows 
a uniform treatment of states with even an odd soliton number. The derivation also 
provides a simple link between the counting function of the Destri-deVega equation 
and the Y-system component functions of the TBA approach. 

After considering the continuum limit of the lattice model, the only input we 
need is the distribution of zeroes of the TBA Y-functions (and their sign at infinity). 
This information can be read off from the infinite volume solution of the problem, 
which we found explicitly. Our main assumption in this paper is that what we find 
for the infinite volume solution remains qualitatively valid for all physical volumes. 
We verified this assumption numerically by comparing the TBA results with those 
obtained by the alternative Destri-dcVcga nonlinear integral equation. 

The simple pattern of the excited state TBA systems found here gives us some 
hope that it will be possible to find TBA equations for excited states also in other 
models, even if no DdV type alternative is available there. 
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A Appendix 

In this appendix we recall the solution \21\ |22] of the basic functional equation 

/(e+o/(e-^) = 5(o. (165) 

Here /(^) is assumed to be meromorphic in the strip |Im,^| < 1 + e, satisfying the 
reality condition 

no = fin- (166) 

Let {ra}a=i be the set of zeroes of /(^) with |ImrQ,| < 1 (not necessarily different 
from each other) and similarly {p/jj^^i the set of (not necessarily different) poles 
with |Imp^| < 1. Moreover, for large (real, positive) ^ 

/(0~^/oee^°«e-^^°^^5« (167) 

asymptotically. Here a = ±1 is a sign, /o > 0, Cq, Ci and A are real. 

The function B{C,) is real analytic in the strip |Im^| < e and non- negative for 
real C,. Its asymptotics is given by 

BiO-^Boe^e'"'^, (168) 

where i?o > and bo, bi are real. Comparing ()167|) with ()168|1 gives 

/o = ci = ^, Co = ^. (169) 

The key observation that allows the transformation of the Y-system functional 
relations, which are of the form p65|) . to TEA integral equations is the fact that 
the solution of ()165p . satisfying all the above requirements is uniquely given by 

11/3=1 - pp) 

where 

r(0 = tanh {^i) (171) 

and 

4 coshfK-t.) ^ ' 

j3{0 is clearly analytic in the strip |Im,^| < 1 and for real u (when B{u) > 0) 

1 ?■ \nB(v) 

I3(u±{) = -\nB(u)±-V dv , ^} ' (173) 
^ ^2 ^ ^ 4 sinhf(t;-M) ^ ^ 

Here V indicates principal value integration. 
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